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1 Introduction 

Technological developments in recent years have 
made it possible to simultaneously record multi¬ 
channel MEG signals with whole-head magneto¬ 
meters, and the relevant research is providing more 
and more impressive results which are enriching our 
understanding of the functions and pathologies of 
the human brain (e.g. [1~5]). One problem in MEG 
signal analysis is that the magnetoencephalographic 
signal is often contaminated by cardiac artifacts 
which can be several times stronger in magnitude 
than spontaneous brain activities [6] (Fig.l). 
Although coherent averaging can reduce the effect 
of cardiac artifacts on event related fields (ERFs), 
the potential effect of such strong artifacts is still a 
concern of many studies based on MEG, particularly 
when studying spontaneous brain activities. 



Fig. 1 One example of MEG recording with cardiac 
contamination 


In this paper, we address this problem by applying 
the common spatial subspace decomposition 
(CSSD) algorithm [7], This method is based on the 
dissociation of topographical subspaces common 
and specific to the MEG measurement matrix and 
the pure artifact matrix. The artifact signal was 
obtained by averaging MEG measurement time- 
locked to the R wave of the cardiac activity. 
Applying CSSD to these two signal matrices can 
decompose the original MEG measurement into two 
parts: one part with the same spatial subspace as that 
of heart beat artifacts, the other part relating to brain 
activities. The artifact can then be removed as the 
common components. 


2 Methods 

2.1 CSSD 

By simultaneous diagonalization of two multi¬ 
channel epochs, CSSD can find common spatial 
factors in both matrices and then decompose the 
components of common spatial factors from the 
original epochs. In the current study, the pure 
artifact epoch (control epoch), containing one 
complete heart beat, was obtained by averaging the 
MEG recording locked to the R peak. The MEG 
recording epoch of arbitrary length was taken as the 
target epoch. Because the cardiac artifact in the 
target epoch and the pure artifact epoch can be 
reasonably treated as having same spatial 
topographies, applying CSSD to these two matrices 
would decompose the artifact from the target epoch. 
As discussed in [7], in order to get perfect 
decomposition by CSSD, the cardiac artifact and the 
brain signal in the MEG recording need to be un¬ 
correlated. This condition can be thought of as 
holding quite well due to the different mechanism of 
electrical current generation in the brain and in the 
heart. 

2.2 Simulation 

Three seconds of signal from a 148-channel (BTi) 
MEG recording from one subject with no obvious 
cardiac artifact was taken as the artifact-free MEG 
signal. One averaged cardiac artifact epoch of 
550 ms from another subject was added to that 
artifact free signal to create simulated target epochs. 
By adding different numbers of heart beats and 
modulating the amplitude of the artifact, 6 different 
target epochs of 3 seconds were created. The time 
points where the artifacts were added and the 
amplitude modulating factors are listed in Table 1. 
CSSD was then applied to the artifact matrix and 
each one of the 6 epochs to remove the artifact from 
the simulated target epochs. In CSSD procedure, 25 
factors were used to account for MEG signal, and 10 
factors for artifacts. 










Table 1 Mixing cardiac artifact with artifact-free 
MEG signal 


epoch 

heart beat number 

modulating factor 

1 

5 at 

149, 727, 1293, 
1932, 2526 ms 

1,0.9, 0.96, 0.8, 0.92 

2 

4 at 

206, 973, 1802, 
2363 ms 

0.9, 0.8, 0.85, 0.88 

3 

3 at 

401, 1291,2145 

ms 

0.85, 0.82, 0.78 

4 

3 at 

637, 1474, 2363 

ms 

0.7, 0.66, 0.69 

5 

2 at 

332, 1618 ms 

0.90, 0.60 

6 

2 at 

637, 1913 ms 

0.6, 0.55 


2.3 Test with experimental data 

MEG recordings from 6 healthy subjects (3 males, 
mean age: 23) taking part in a language experiment 
were analyzed. These recordings had obvious 
cardiac artifact contamination. MEG signal was 
recorded at a sample rate of 508.63 Hz with online 
filtering of 0.1-100 Hz. 

First, triggers were set to the R peak of each artifact 
complex by band-pass filtering (5-20 Hz) and 
maximum detection. The pure artifact was then 
obtained by averaging artifact epochs (-150 ms to 
400 ms with respect to triggers). Prior to averaging, 
epochs were screened for eye movements and other 
artifacts. The epochs were excluded from averaging 
if (a) the standard deviation of the EOG within a 
sliding 200-ms time window exceeded 30 pV, or (b) 
the standard deviation of any MEG recording 
channel within a sliding 3-second time window 
exceeded 1100 fT, or (c) the drift of the MEG 
recording within the epoch was over 3000 IT. 
Second, the MEG recording between two 
consecutive rejected epochs was taken as the target 
epoch. CSSD was then applied to the target epoch 
and the pure artifact epoch to remove the cardiac 
contamination from the MEG recording. The result 
is the ‘artifact-free’ MEG recording. 

Finally, the resultant ‘artifact-free’ MEG recording 
locked to the former R peak of heart beat artifacts 
was averaged to get the reduced pure artifact. 
Ideally, the resultant artifact signal would be 
negligible compared to the original pure artifact 
signal. 


2.4 Source localization with a phantom 

To test the effects of the current algorithm on source 
localization, a phantom experiment was conducted. 
As a phantom head, a saline-filled sphere provided 
by the manufacturer of the MEG system was used. 
Five artificial current dipoles at different depths 
within the phantom were used to generate artificial 
MEG signals. The radius of the phantom sphere was 
7.61 cm. The locations of the dipoles within the 
phantom with respect to the ‘nasion-ear’ coordinate 
system are listed in Table 2. The dipoles were 
activated in turn and the current through the dipoles 
was kept unchanged for different dipoles. For each 
source, 500 epochs were created, consisting of 
200 ms signal free data and 200 ms of 10 Hz sine 
wave. 


Table 2 Parameters of the planted dipoles 


dipole 

No. 

location (cm) 

length 

(cm) 

X 

y 

z 

1 

-0.17 

1.69 

6.05 

0.539 

2 

-0.13 

1.39 

5.09 

0.517 

3 

-0.11 

1.12 

4.13 

0.531 

4 

-0.09 

0.86 

3.16 

0.546 

5 

-0.07 

0.61 

2.2 

0.533 


To simulate the cardiac contamination, MEG 
recording from one subject with obvious heart beat 
artifacts was mixed with the MEG recording from 
the phantom. The amplitude of the MEG signal from 
the subject was adjusted to get different signal-to- 
noise ratios. The procedure described in section 2.3 
was used to remove the artifact from the mixed 
MEG recording. For each artificial dipole, the 
averaged phantom signals (from 400 epochs) before 
and after artifact reduction were then fed into 
Curry® (NEURO Scan Labs) for dipole fitting. For 
the dipole modeling, MEG signals were co¬ 
registered with a sphere model based on nasion-ear 
landnarks. 

3 Results 

3.1 Simulation 

The simulation results were shown in Fig. 2. The 
spatial correlation coefficient (CF) is that between 
two topographies at one time point. The temporal 
CF is that between two time waveforms of one same 
channel. The spatial CF depicted in this figure is the 
mean value across the whole time window, and the 
temporal CF in the figure is the mean value across 
all the 148 channels. It can be clearly seen that after 





removing the artifacts, all the CFs are above 0.96. 
The distortion caused by artifacts was effectively 
removed. 
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Fig. 2 Improvement of correlation coefficients 
between the artifact-free MEG signal and the 
artifact contaminated MEG before (blue) and after 
(red) CSSD. The horizontal axes refer to the 6 cases 
in Table 1. Left: spatial aspect; Right: temporal 
aspect. 


3,2 Experimental data analysis 



Fig. 3 Pure artifact signals before (blue) and after 
(red) CSSD. About 600 trials were averaged for 
each subject. 


Fig.3 shows one channel of the cardiac artifact 
before and after artifact reduction by CSSD for each 


of 6 subjects. Compared to the original artifact 
signal (blue), the averaged artifact after reduction 
(red) is negligible. The reduction ratios between the 
peak values of reduced and non-reduced artifacts 
are: (left column, top town) 1/40, 1/52, 1/26; (right 
column, top down) 1/63, 1/25, 1/91. The mean 
reduction ratio for these 6 subjects is 1/50. The 
reduction is very effective. 

3.3 Localization with the phantom dipoles 

The dipole localization results with the phantom 
experiment are depicted in Fig. 4. The signal-to- 
noise ratio (SNR) was estimated in Curry using the 
200 ms signal free interval in the phantom 
experiment as the noise time window. The SNR 
value shown in this figure is that without artifact 
reduction. In the present study, after the artifact 
reduction by CSSD, SNR was improved by a factor 
of between 2 and 5. When localizing the dipoles in 
Curry, the regularization procedure was applied to 
compensate the effect of noise. 

The localization accuracy improvement was defined 
as: mislocation before CSSD - mislocation after 
CSSD. As shown in Fig. 4, the localization accuracy 
can be improved by 3-4 mm for superficial dipoles 
(No. 1&2), 5-6 mm for sources at moderate depths 
(No. 3&4), 6-9 mm for the deep source (No. 5). 



SNR 


Fig. 4. Dipole localization accuracy improvement 
by the heart beat artifact reduction. 

4 Discussion 

Through simulation, experimental data analysis, and 
the phantom dipole localization, we have 
demonstrated the applicability and effectiveness of 
CSSD in removing the cardiac artifact from MEG 
recordings. 

































Because CSSD is based on dissociation of spatial 
subspaces, the target and control epochs in current 
study can have different time lengths, and the 
artifacts in different epochs can have totally 
different magnitude. As long as the artifacts in these 
two epochs share the same topographic subspace, 
they can be removed by the present algorithm. 

One concern over the present algorithm is the 
possible distortion to the topographies of brain 
activities. As discussed in [7], such a distortion in 
CSSD is caused by the correlation between artifact 
and brain activities. Since a high correlation 
between the heart beat artifact and the brain activity 
is not anticipated, the distortion could be negligible. 
At least, the simulation result in Fig 2 demonstrated 
that the spatial distortion caused by artifacts can be 
effectively reduced by the present method. 

The cardiac artifact in MEG is generated by a very 
deep source. Therefore, the effect of the artifact on 
the localization of deep sources will be more 
significant than on the superficial sources. That is 
why the localization accuracy improvement for the 
deep source in the phantom signal analysis is larger 
than for others. The improvement of the localization 
accuracy also demonstrates that possible distortions 
to the topographies of interest caused by CSSD are 
negligible. 

In the phantom experiment, the artificial dipoles in 
the phantom were located at about the center point 
with respect to the left and right axis. Since MEG 
sensors at the edge of the dewer usually pick up the 
strongest artifact signal, the effect of the artifact on 
source localization may be more serious when the 
source is located to the more left or to the more right 
(e.g. when temporal brain regions become 
activated). 
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